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LL£      (  tX  ABSTRACT 


A  study  is  made  of  the  circumstances  under  which  there  will 
be  incipient  buckling  in  an  infinite  plate  containing  a  circular 
hole  and  subjected  to  uniform  uniaxial  tension  at  a  great  distance 
from  the  hole.   This  buckling  is  a  consequence  of  the  fact  that 
near  the  hole  and  in  line  with  the  direction  of  the  applied  tensile 
loading,  there  are  areas  in  which  both  principal  membrane  stresses 
are  compressive.   Values  of  critical  tensile^  stress  are  given, 
together  with  expressions  for  the  corresponding  lateral  deflection 
patterns.   The  results  are  in  the  form 

S   =  KE(t/a)2 
cr 

where  E  is  Young's  modulus  of  elasticity,  t  is  plate  thickness,  a 
is  hole  radius,  and  K  is  a  buckling  coefficient  having  the  following 
values  for  the  first  four  modes:   7.65,  7.80,  11.7,  12.3,  A 
comparison  is  made  with  some  experimental  determinations  of  Danis 
for  plates  of  finite  width. 
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y*  A  compacted  notation  for  Y„,  in  Table  5  with  I  and  J 

NM 
replacing  N  and  M  respectively 

y  A  compacted  notation  for  Y   in  Table  5  with  K  and  L 

replacing  N  and  M  respectively;  also  a  coordinate  of 
the  Cartesian  system;  also  a  scalar  parameter  in  an 
eigenvalue  problem  in  Appendix  D 

An  unknown  coefficient  which  is  used  to  satisfy  the 
boundary  conditions  for  each  term  in  N  or  M 

The  number  of  terms  taken  in  N  in  the  series  for  W 

The  number  of  terms  taken  in  M  in  the  series  for  W 

A  scalar  parameter  in  eigenvalue  problems 

Normal  stress 

Shear  stress 

Q  An  angular  coordinate  in  the  jbylindrical  system 

-,)  Poisson's  ratio  (v  =  .3  is  used  throughout) 

OTHER  NOTATION 

A  double  underline  indicates  a  square  matrix.  (A) 

A  single  underline  indicates  a  column  vector.  (A) 


11. 


Y 

NM 

(X 

P 

A 

• 

^°3l°"r, 

^ 

X©, 

T 

ACKNOWLEDGEMENTS 

The  author  takes  this  opportunity  to  thank  Professor  John  E.  Brock, 
for  his  patience,  understanding,  and  guidance  in  this  endeavor;  and  to 
thank  Professor  Robert  E.  Newton  and  Professor  Gilles  Cantin  for  their 
helpful  suggestions. 


12. 


1.   INTRODUCTION 

The  problem  of  determining  the  stress  situation  near  a  hole 
in  a  plate  under  tension  (see  Figure  1)  has  a  well  known  solution 
which  is  of  great  practical  importance. 

The  original  work  on  this  problem  was  done  by  Kirsch  (19), 

who  solved  the  problem  of  the  stress  distribution  in  an  infinite 

plate  pierced  with  a  round  hole  and  subject  to  uniaxial  tension 

i 

at  a  great  distance  from  the  hole.   Howland  (17)  treated  a  similar 
problem  with  an  infinite  strip  of  finite  width.   Savin  (5)  and 
Hogan  (15)  have  compiled  solutions  produced  by  other  authors  for 
variations  in  hole  geometry,  as  square,  triangular,  etc.,  for  both 
finite  and  infinite  plates.   Lur'e  (3)  has  solved  the  problem  with 
a  hole  in  a  shell  of  uniform  curvature. 

An  important  result  obtained  from  these  solutions  is  that  the 
stresses  near  the  hole  become  much  greater  than  those  imposed  upon 
the  plate  at  a  distance  from  the  hole.   This  has  given  rise  to  the 
wide  publication  of  many  stress  concentration  factors  which  are 
important  in  design  work  and  which  relate  the  maximum  stress  to  the 
nominal  stress  at  a  distance  from  the  hole. 

Another  feature  of  the  stress  situation,  which  is  not  so  widely 
recognized,  is  that  certain  areas  of  the  plate  are  under  uniaxial  or 
biaxial  compression,  although  tension  is  imposed  externally.  These 
areas  of  compression,  given  the  proper  conditions  of  plate  thickness 


*  These  terms  will  be  employed  herein  to  mean  that  one  or  both 
principal  stresses  are  compressive.   Specifically,  "uniaxial  compression" 
means  that  one  principal  stress  is  compressive  and  one  is  tensile; 
"biaxial  compression"  means  that  both  are  compressive. 
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and  external  tensile  load,  will  buckle,  possibly  causing  failure  of 
the  plate.   Houghton  and  Rothwell  (16)  have  shown  that  in  plates 
with  holes  having  a  very  small  gross  curvature   the  stress  concen- 
tration factor  can  increase  appreciably.   It  is  reasonable  to  assume 
then,  that  even  though  local  distortions  of  a  plate  may  be  acceptable, 
such  distortions  could  increase  the  stress  concentration  beyond  that 
expected  for  a  flat  plate.   Failure  of  the  plate  by  yielding  could 
then  occur. 

There  is  no  known  investigation  of  this  buckling  problem,  other 
than  an  experimental  investigation  by  Danis  (14),  to  determine  the 
critical  stress  for  the  buckling  of  a  finite  plate  with  a  round  hole. 

This  paper  presents  a  theoretical  investigation  of  the  incipient 
buckling  problem  for  an  infinite  plate  containing  a  round  hole.   This 
problem  is  of  interest  in  itself  and  for  the  practical  reasons  given 
in  the  preceding  discussion.   Results  are  presented  on  page  33. 
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2.   PROBLEM  DESCRIPTION 

Consider  a  plate  infinitely  extended  in  two  directions  and 
pierced  with  a  hole,  as  shown  in  Figure  1. 

s 

M  U  M  t  ♦  n  M  I   t   > 


Figure  1 
x  and  y  are  Cartesian  coordinates,  r  and  0,  cylindrical  coordinates. 
The  z  axis,  not  shown,  is  perpendicular  to  the  paper  at  the  center 
of  the  hole,  and  extends  in  the  positive  direction  toward  the  reader. 
Consider  all  quantities  to  be  in  units  of  hole  radius,  a,  which 
essentially  reduces  them  to  dimensionless  quantities,  with  hole  radius 
equal  to  one. 

The  plate  is  placed  in  uniform  tension,  S  psi,  at  infinity  in 
the  x  direction.   The  solution  of  the  stress  distribution  problem  for 
this  plate,  as  originally  solved  by  Kirsch  (19),  and  clearly  presented 
by  Timoshenko  (7)  is: 
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crr  = 

S  ' 

2 

Ce  = 

S 

2 

X9  = 

S 
2 

/-p>)  t  (/+pf-fzl 


1-^  +  ^  S//V2S 


(1) 
(2) 
(3) 


This  set  of  equations  is  henceforth  referred  to  as  the  Kirsch 
solution. 

It  is  observed  that  the  maximum  tensile  stress  is  3  times  the 
stress  applied  at  infinity,  and  occurs  at  r  =  1,  and  9  =  +7T/2. 
The  maximum  compressive  stress,  equal  in  magnitude  to  the  applied 
tensile  stress,  occurs  at  r  =  1,  and  9  =  +7f.      From  this  it  can 
be  seen  that  if  buckling  were  to  occur  in  the  plate,  it  would  be 
most  apparent  near  the  hole  along  the  x  axis. 

A  better  understanding  of  the  buckling  problem  is  obtained  by 
considering  the  minimum  principal  stresses  (most  negative,  largest 
compressive).   From  the  Kirsch  solution,  the  minimum  principal 
stresses  are  found  and  shown  in  Figures  2  and  3. 

Figure  2  shows  one  quarter  of  the  plate,  with  the  hole  centered 
at  the  top  left.   The  stress  distribution  is  symmetrical  with  respect 
to  the  x  and  y  axes.   This  picture  was  prepared  with  consideration 
only  to  whether  the  stress  situation  at  a  particular  point  is  biaxial 
tensile,  biaxial  compressive,  or  uniaxial  compressive. 

Figure  3  shows  the  same  section  of  plate  as  Figure  2.   This 
presentation  maps  the  relative  values  of  the  algebraic  minimum 
principal  stress  at  each  point. 
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FIG  3  -    MAP    OF   THE    MINIMUM   PRINCIF&L   STRESS 
AROUND  A  HOLE 
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It  can  be  seen  from  these  figures  that  there  is  an  area  having 
biaxial  compression.   Accordingly,  buckling  is  expected  to  take  place 
for  some  sufficiently  great  value  of  applied  load. 

Numerous  methods  have  been  developed  for  solutions  of  plate 
buckling  problems.   In  general,  these  methods  involve  either  the 
satisfaction  of  the  equilibrium  equations  or  minimizing  the  strain 
energy  of  the  plate  (Energy  Methods).   A  compilation  of  many  solutions 
to  plate  buckling  problems  was  made  by  Gerard  and  Becker  (11).   This 
is  an  excellent  reference  for  both  method  and  practical  solutions. 

One  of  the  first  energy  methods  was  developed  by  Ritz  (20),  and 
independently,  by  Lord  Rayleigh  (4).   This  method  is  presented  in  a 
rigorous  manner  by  Sokolnikoff  (6).   Budiansky  and  Hu  (12)  give  a 
detailed  example  of  the  application  of  this  method  to  a  plate  buckling 
problem. 

The  Rayleigh-Ritz  method  involves  determining  the  total  strain 
energy  in  the  plate  as  a  function  of  an  assumed  buckled  configuration. 
Equilibrium,  and  the  stable  buckled  deflection,  W(x,y),  exists  when 
this  strain  energy,  U  ,  is  minimized  with  respect  to  W(x,y). 

To  do  this,  first,  a  family  of  functions,  Q   ,  are  assumed,  which 
are  dependent  upon  parameters  A  ,  such  that  the  buckled  deflection, 
W(x,y),  can  be  represented  by: 

W(*,»)  =  fK^A,,^,  •••)  (4) 

where  ^  is  such  that  the  boundary  conditions  on  the  plate  are 

satisfied  for  any  values  of  the  A. . 

l 

The  functions  S>  axe   substituted  into  the  strain  energy  expression 

U   so  that 
t 

Ut=jfjf#t.  §(*,y>\bX)'  •  0    tedy  (5) 

where  U   is  the  integrated  total  energy  in  the  plate. 
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In  order  to  minimize  the  energy  function,  the  following  set 


of  equations  is  formed 


df\, 


o 


-  o 


(6) 


If  this  set  of  equations  has  solutions  A  A   . ..,  the  function 

represents  the  buckled  equilibrium  condition,  the  boundary 
conditions  of  the  problem  being  satisfied  by  ^  for  any  values  of  the 

A.. 

1 
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3.   BOUNDARY  CONDITIONS 

For  the  problem  of  an  infinite  plate  with  a  hole,  two  boundaries 
are  of  concern.   At  infinity,  deflection,  slope,  and  curvature  are 
zero.   Around  the  hole,  a  free  edge  condition  exists. 

The  Kirsch  solution  is  used,  which  already  implies  zero  tractions 
around  the  hole,  and  uniaxial  tension  at  infinity,  so  these  conditions 
have  already  been  satisfied. 

For  a  free  edge  to  exist  around  the  hole  means  that  there  is  no 
shear,  twist,  or  moment  applied  at  the  edge  of  the  hole.   Following 
Timoshenko  (8),  and  considering  the  cylindrical  coordinate  system 
established  in  Figure  1,  the  conditions  for  a  free  edge  are 

f  Q  \  =0  (shear  =  0)  (7) 

(  MYt)         =0   '  (twist  =  0)  (8) 

'Y  —  \ 

(Mr)^.   =  O  (moment  =  0)  (9) 

Kirchoff  (18)  has  shown  that  the  twist  condition  can  be  reduced 
to  a  statically  equivalent  shear  force.   The  original  requirements 
of  zero  shear  and  zero  twist  can  be  reduced  to  one  requirement  such 
that  the  sum  of  the  shear  forces  and  the  static  equivalent  of  the 
twisting  moments  equals  zero.   The  conditions  then  become 


(Qv- 


3Mrt 


=  o 


r-\ 


►H.,-0 


(Kirchoff  shear 
(moment  =0) 


=  0) 


where  Q    ,   M   ,    and  M        are   expressed  as 
xr        r  rt 


(10) 
(ID 


Qr=  -D& 


Mr  =  -D 


Mvt=  W'-°>  T 


aw       I     3W       i     dxW 
Br1         r  '  d\r  +     rx'  a#>- 


+  -o[±.|w.  J.  £w 

1 

-■o) 

i     92w        i      aw 

(12) 
(13) 
(14) 
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When  the  appropriate  differentiations  and  substitions  are  made, 
the  boundary  condition  equations  become 


Jw._L.gtW   j_  3W   (2--0)  9s  W   (3-t))  3aW 
r*   «r  9y-2   r2*9r  "^  r3-  'arS©3--  v-3  'a©2 


3 


(Kirchoff 

Shear  =  0) 


=  0 
V  =  i 


9r2    r  1  dr   ^'d©2 


=  0 


(Moment  =  0) 


=  O 
y-=oo 


(15) 
(16) 

(17) 


As  previously  noted,  each  member  of  the  family  of  functions 
must  satisfy  the  boundary  conditions  in  the  Ray leigh-Ritz  method. 
For  a  choice  of  ^  let 


oo  oo 

r-1  v^  -n 


N  =  l  M=o 


MM 


ANMcos  Me 


0\r 


Bnmsin  MS 


(18) 


where  A   and  B._,  are  arbitrary  coefficients  that  are  varied  to 

NM      NM 

minimize  strain  energy  in  the  plate. 

The  boundary  conditions  of  eq.  17  are  satisfied  by  this  series 

since  at  r  =  o°^W  and  all  of  its  derivatives  are  zero.   Conditions  of 

equations  15  and  16  are  satisfied  by  substituting  equation  18  into 

these  two  boundary  condition  equations  (see  Appendix  A).   Such  a 

substitution  results  in  two  equations,  15  and  16,  with  two  unknowns, 

XXT„  AND  Y._„,  for  each  conbination  of  N  and  M.  A  unique  value  of 
NM      NM  M 

* 

X  „,  and  Y„_,  is  found  which  allows  each  term  of  equation  18  to 

NM      NM  n 

satisfy  the  free  edge  boundary  conditions  for  any  N  and  M  chosen. 
From  Figures  1  and  3,  it  is  seen  that  there  is  symmetry  with 
respect  to  both  the  x  and  y  axes.   This  means  that  when  W  is  written 
as  a  Fourier  series  in  0  (eq.  18),  the  series  contains  either  cosine 


*  A  term  in  this  sense  means  the  full  expression 

ANMcos  M0 


-n 


+  Xnm+  Ynm  . 

v    rx 


ov- 

6mm  sin  Me 
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terms  or  sine  terms,  but  not  both.   Additionally,  the  requirement 
of  symmetry  around  both  the  x  and  y  axes  means  that  only  odd  or 
even  terms  in  M  need  be  taken.   Thus  the  Fourier  portion  of 
equation  18  reduces  to: 

A_  cos  M0,  with  M  =  1,3,5,... 

or 

A„,  cos  M0,  with  M  =  0,2,4,... 
NM 

B._,  sin  M9,  with  M  =  1,3,5,... 

NM 

B   sin  M9,  with  M  =  2,4,6,... 

This  gives  four  separate  expressions  for  W,  one  for  each 
possible  condition  of  symmetry.   The  full  Fourier  series  could  be 
taken,  but  if  this  is  done,  three  quarters  of  the  coefficients,  A^ 
or  B   are  zero  and  considerable  computer  storage  space  is  wasted. 
It  is  shown  later  than  each  of  these  series  represents  a  different 
buckling  mode. 


or 


or 
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4.   ENERGY  IN  THE  PLATE 


To  this  point,  a.  family  of  functions  W(r,9)  (equation  18)   have 
been  developed,  which  satisfy  the  boundary  conditions  independently 

of  the  coefficients,  A„,  or  B._„,  and  are  adequate  to  represent 

NM     NM  n  r 

possible  shapes  of  the  plate.   Attention  is  now  turned  upon  formu- 
lating an  expression  for  the  total  strain  energy,  and  then  minimizing 
this  expression  to  obtain  the  critical  buckling  stress. 

The  energy  in  the  plate  can  be  divided  into  two  forms.   First, 
there  is  the  membrane  energy  which  is  the  energy  stored  because  the 
mid-plane  is  stretched.   Second,  there  is  the  bending  and  twisting 
energy  which  results  from  the  plate  buckling  out  of  its  plane. 

A  basic  assumption  is  that  the  plate  deflection,  W(x,y),  is 
inf initesimally  small,  and  that  the  mid-plane  strains  may  be  re- 
presented by 


e*~   a*+  2  1571 
_  3y     ±  (awf 

e3"   5J+    2   lay   ) 

r  -    ^±+  9v  f  aw  aw 

r*3~     Sy      dx  r  d%  '  3g 


\ 


(19) 


/ 


Terms  of  higher  order  in  W  than  these  are  neglected,  since  W 
is  infinitesimal;  those  terms  which  are  retained  are  essential  for 
the  calculation  of  bending  energy. 

In  the  above  expressions  £~}  etc.,  represent  the  stretching 
of  the  plate  in  its  own  plane.   The  terms  _L (  —^ j ^    etc., 
represent  the  difference  between  a  length  of  flat  plate,  and  the 
corresponding  straight  line  length  when  the  plate  is  bent  (difference 
between  arc  and  chord) . 
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Stresses  in  the  x  and  y  directions  are  composed  of 
(a)   Membrane  stresses  of  finite  magnitude,  which  come 
from  the  Kirsch  Solution  and  are  denoted  by  C^  Q"L  and 
Tx         below;  and   (b)   Infinitesimal  bending  and  twisting 
stresses  which  do  not  appear  explicitly  below,  but  which 
are  involved  in  the  bending  and  twisting  energy. 
The  normal  forces  per  unit  length  (N  ,N  )  can  be  determined  from 
N  =  0"^n  ,  etc.   The  membrane  energy  in  the  plate  is 


MEM 


\  J/(Nxex  +  N3ey  +  N^r^jckcty 


(20) 


Substituting  the  strains  (eqns.  19),  and  letting 


UMEM  "  Ul  +  U2 


where 


<ixty 


(21) 


(22) 


and 


a-  fsw\2      lawi2,  r  .  aw.  awl  jv  ju 


(23) 


an  expression  for  U,_.,  is  found  in  terms  of  the  stresses,  strains, 

r  MEM 

and  the  deflection. 

It  can  be  seen  that  the  membrane  energy  is  composed  of  two 

terms.   One,  IL  ,  relates  the  work  done  by  the  stresses,  O"-^  ,  etc, 

through  a  distance  2—   ,  etc.   The  other  term,  U_ ,  relates  the 

ay.  *- 

I  /  9W\2 
work  done  by  the  stresses,  0"x  ,  etc.,  through  a  distance  r-(  ^rrj) 

etc.,  as  the  plate  bends  out  of  plane. 

At  this  point,  the  assumption  is  made  that  the  Kirsch 

solution  is  still  valid,  even  if  the  plate  has  buckled  infinitesimally . 

This  means  that  U   is  independent  of  the  lateral  deflection,  W,  and  can 

be  treated  as  a  constant.   This  assumption  is  justified  by  considering 
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that  the  problem  is  one  of  incipient  buckling,  not  postbuckling. 
This  assumption  considerably  simplifies  the  problem,  because  it 
uncouples  the  membrane  stress  distribution  from  the  buckled 
deflection. 

An  expression  for  the  bending  and  twisting  energy,  U„,  is 
obtained  by  considering  the  moments  in  the  plate  and  the  respective 
rotations  at  a  point  as  functions  of  plate  deflection.   Timoshenko  (9) 
shows  this  expression  to  be; 


dx3-     9<j* 


-U-*) 


alw  b2w  \dz\N 


9X1  By3-   l3x9y 


dx<J3 


(24) 


Adopting  the  notation  Wy-=  =r—  >  Wg=  ^^ \   etc.^   and  transforming 
the  previous  energy  equations  to  cylindrical  coordinates,  with  the 
proper  limits  of  integration  for  the  problem  of  a  hole  in  an  infinite 
plate,  equations  22,  23,  and  24  become,  respectively 


U   =  CONSTANT 
o       I 


[vrVJ?+^  +  lrt^kV*L  )^rd9 


(25) 
(26) 


»3-f 


-277V  °° 


o  -,\ 


W^+^  +  ^ 


r* 


-(2-0) 


w^|V^  +  V^e|        (27) 


v3-   r 


Hrd0 


CTy-^e  and  'TyQ  are   obtained  from  the  Kirsch  solution  (eqns.  1,2,  and  3). 
W  and  its  derivatives  are  obtained  from  the  assumed  series  (eq  18). 
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5.   SOLUTION  OF  THE  PROBLEM 

With  the  expression  for  total  energy  (U  =  U  +  U_  +  U  ) 
formulated,  the  next  step  is  to  substitute  the  Kirsch  solution, 
and  the  assumed  series  for  deflection  and  its  derivatives.   Then 
this  expression  must  be  minimized  with  respect  to  each  of  the 
unknown  coefficients,  A^M  or  B   .   The  substitutions  and  inte- 
grations are  shown  in  Appendices  B,  C,  and  D.   It  is  sufficient 
for  our  purpose  here  to  say  that  what  arises  from  the  integration 
of  the  expressions  for  U9  and  U  ,  are  two  quadruple  series,  with 
indices  I  and  K  summed  over  the  range  of  N,  and  J  and  L  summed 
over  the  range  of  M.   The  quadruple  series  stems  from  the  multi- 
plication  of  two  double  sums.   Equations  26  and  27  now  become 

U2  =   f££EX   GWA„-AKL  (28) 

U3=    irLEII    Px«,"Aw-AKt  (29) 


2 

I   J   K  U 


Eh3 


where  the  flexural  rigidity  of  the  plate,  D,  equals  — ; -=t- . 

IZ(/--02) 

Q     and  P     are  different  for  each  value  of  I,  J,  K,  or  L. 

1JK-L        1JK.L 

S  is  the  stress  applied  to  the  plate  at  infinity   and  is  factored 

from  the  Kirsch  solution  after  those  expressions  are  substituted 

into  the  energy  integrals.   A   and  A   are  the  unknown  coefficients 

IJ       KL 

in  the  series  expression  for  deflection  (eq.  18). 

Following  the  Ray leigh-Ritz  procedure,  U   is  differentiated 
with  respect  to  each  of  the  unknown  coefficients,  so  that  the 


The  unknown  coefficients  for  the  cosine  series,  A^,  only  will  be 
used  henceforth,  although  these  could  be  the  sine  series  coefficients, 

B„,,  equally  as  well. 

NM 
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following  set   of   simultaneous   equations   result   i        \ 


3A12 


+ 


§^„„  +  £A 


2  n/ii 


A„ 

A 


(30) 


12 


/ 


The  last  term  in  each  equation  arises  from  the  fact  that 

As  a  practical  matter,  the  equations  arrived  at  in  this 

investigation  are  too  complicated  to  permit  considering  the 

use  of  a  doubly  infinite  set  of  terms  in  the  double  series. 

Accordingly,  it  is  henceforth  presumed  that  truncated  series 

are  used  and  that  there  is  but  a  finite  number  of  unknown  A„„ ,'s. 

NM 

Simplification  in  notation  is  made  at  this  point  if  the 
quadruple  subscripts  are  replaced  by  a  double  subscript  arrange- 
ment. 

Let  B  equal  the  number  of  terms  taken  in  M,  and  o(   equal 
the  number  of  terms  taken  in  N.   Define 

i'=(i-0p  +J 

J'=(K-I)p  +  L 

Using  these  equations  a  unique  pair  of  subscripts  can  be  changed 

to  a  unique  single  subscript,  so  that  A   =  A  ,,  and  Q     =  QT  ,  , , 

1J     i  1JKL     i  J 

etc. 

For  example,  if  4  terms  were  taken  in  M,  and  we  wished  to 

convert  Q^/no  to  some  doubly  subscripted  Q  ,  ,,  then  1=7,  J  =  4, 
7  4y  Z  JL  J 


K  =  9,  and  L  =  2;  and 

I'  =  (7-1)  4+4, 

=  28 

so  that  Q-,,00  becomes  Q 
/  4-yz 


(28) (34) 


J1  =  (9-1)  4  +  2 

=  34 
Similarly 
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A, i            A, 

Qmi   -* 

Qn 

An.  ""*'    Aj. 

Qia.\i    — ► 

Qai 

An  -  A3 

61  IM2          * 

Q.a 

etc. 

etc. 

Also, 

let 

J=  (*-l)'P 

+  ? 

-' 

o(. 

P 

where 

J  is    the    large 

.St 

value 

of 

I1 

and  J'    respectively. 

Substituting 


D  =  - 


Etf 


and 


9Ut 

8A, 


O  - 


Fk1  Ki3-'    ^     'u' 


9^  =  0  = 
3A,. 


J' 


J' 


|6S(l-0a)\p         ,    T    ' 

I       £j?  JK2T'    +      '2T' 


^A5 


|  6Sfl--0t))n  +T 


'   -J     t' 


i7j 


/2(/--0a) 

RlV  =  QJW      and        Tr'i«  =  PrV  lf 

RiV  =  2-Qr'j'      «*d       Trr  =  2-PxV      If       r'  =  J' 
equation   30   simplifies    to 


\ 


J' 


(31) 


When  the  substitution 

\  =   fcS(i-T)1)  (32) 

A     Eh2 
is  made,  and  the  above  set  of  equations  expressed  in  matrix 

notation,  the  following  matrix  eigenvalue  problem  results 

(RX  +  I)A=0  (33) 

where 
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R  = 


"21        Ril 


R31    Rya 


R23 


R55 


(34) 


T  = 


T11     T,z 

T2l      "Htt 


Tj,    Tja- 


T,ff 

T2J 


T55 


(35) 


and 


A  = 


r      n 
A, 
A* 


A5 


(36) 


The  matrix  T  is  positive  definite,  since  it  arose  from  the 
bending  energy  integral.   The  R  matrix  arose  from  the  membrane 
energy  integral,  which  contains  terms  from  the  Kirsch  solution. 
This  matrix  is  not  necessarily  positive  definite.   Both  matrices 
are  symmetric.   These  properties  are  important  to  the  method 
utilized  in  solving  the  matrix  eigenvalue  problem  (see  Appendix  D ) . 

The  question  of  whether  differentiating  U  and  setting  these 
derivatives  equal  to  zero,  produces  a  minimum  or  a  maximum,  has  not 
been  discussed.   Either  situation  can  occur  in  the  problem  as  it 
has  been  developed.   Some  of  the  eigenvalues  obtained  in  a  solution 
of  equation  33  represent  maxima  which  are  not  relevant  to  the  plate 
buckling  problem.   The  criterion  for  rejecting  such  answers  is  the 
set  of  equations 
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3H4 
3Af 


(37) 


O 


/ 


If  the  set  of  equations  37  is  satisfied,  then  U   is  minimized. 
When  the  set  of  equations  31  is  again  differentiated  with 
respect  to  the  appropriate  coefficient,  the  set  of  equations  37 
become 


> 


(38) 


daUt  • 


=  XRys+Tjf  =*o 


where  Rtitt  and  ^  ?  i  are  tne  diagonal  elements  of  the  matrices 
R  and  Tj  respectively.  All  of  these  terms  are  positive  numbers, 
since  they  arise  from  the  "squares"  in  the  energy  integrals  (i.e., 
no  cross  product  terms). 

By  rearranging  equation  38  so  that  it  becomes 


*  0"  <  J 


(39) 


It  is  seen  that  for  some  negative  values  of  A  ,  and  all  positive 
values,  the  solution  of  equation  33  represents  a  minimization  of 
energy. 

We  are  interested  in  only  positive  values,  so  that  for  any 
answer  of  interest,  a  minimization  is  assured.   It  should  be  noted 
that  the  negative  values  of  A  that  are  admissible  stem  from 
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buckling  modes  in  the  plate  that  correspond  to  an  applied  compressive 
loading.   These  modes  arise  from  the  reversal  of  sign  in  the  Kirsch 
solution  and  the  associated  reversed  stress  distribution  around  the 
hole.   Modes  of  this  type  should  be  expected  from  a  theoretical 
standpoint. 

When  a  nontrivial  solution  has  been  found  for  equation  33, 
giving  an  eigenvalue  A  and  a  corresponding  eigenvector  A,  then  the 
buckling  stress  may  be  found  by  use  of  equation  32  and  the  buckled 
mode  shape  from  equation  18. 
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6.   RESULTS  AND  CONCLUSIONS 

The  results  of  this  thesis  may  be  presented  compactly  in  the 
form  of  a  formula  for  incipient  buckling,  derived  from  equation  37, 
viz . 

SCf=KE(|)  (40) 

where  S   is  the  value  of  S,  the  tension  at  infinity,  which  causes 
incipient  buckling,  t  is  plate  thickness,  a  is  the  hole  radius  and 
E  is  Young's  modulus  of  elasticity.   Values  of  the  buckling  coef- 
ficient, K,  are  given  in  table  1. 

Values  of  K  were  obtained  from  a  digital  computer  solution  (see 
Appendix  F) ,  for  the  particular  case 

E  =  30  x  10  psi   (Young's  Modulus) 

-u  =  .3  (Poisson's  ratio) 

h  =  .001  hole  radii 

Equation  40  permits  prediction  of  S   for  other  values  of  t/a 
and  E.  However,  the  effect  of  different  values  of  -d   cannot  be  studied 
easily  since  -0  appears  in  the  bending  energy  integral  and  in  the 
boundary  conditions.  All  results  given  are  f or -0  =  .3,  but  the 
computer  program  can  be  modified  by  changing  card  listed  as  0011 
on  page  84 . 

Four  distinct  mode  shapes  are  apparent,  as  determined  by 
different  buckling  coefficients  appearing  for  each.   A  tabulation 
of  K  is  shown  in  Table  1  below. 
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Table 

1 

Values  of 

the 

Buckling 

Coefficient, 

K 

MODE 

K 

1 

7.65 

2 

7.80 

3 

11.7 

4 

12.3 

Since  the  energy  method  involves  the  use  of  a  series  which 
cannot  be  made  infinite  in  the  present  method,  the  actual  buckling 
coefficients  obtained  will  always  be  higher  than  the  true  value. 
This  is  because  the  deflection  of  the  plate  is  only  approximated, 
and  the  energy  never  quite  minimized.  A  solution  with  less  plate 
energy  could  be  found,  if  more  terms  were  taken  in  the  series  for 
deflection. 

Different  values  of  the  buckling  coefficient,  K,  actually 
calculated  by  use  of  the  digital  computer  program  are  shown  in 
Table  2,  and  are  plotted  against  the  number  of  terms  taken  in 
both  the  N  and  M  parameters  of  the  deflection  series,  equation  18. 
These  plots  are  shown  in  Figures  4a  through  4e.   The  projected 
values  listed  in  Table  1  are  taken  graphically  from  these  figures 
as  the  asymptotes  which  the  curves  for  buckling  coefficient  appear 
to  approach,  when  the  number  of  terms  in  N  and  M  increase.   These 
projected  values  are  considered  to  be  the  best  evaluation  of  the 
buckling  coefficient. 

Figure  4f  is  taken  from  the  experimental  data  of  Danis  (14). 
All  of  his  points  for  a  finite  plate  are  plotted  as  K  versus  d/b, 
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Table  2 

Values  of  the  Buckling  Coefficient,  K,  for  Various 
Numbers  of  Terms  in  the  Series 


Mode  1 


6 
8 
10 
12 
14 
16 
18 
20 


a 

•r-l 
CO 

B 

cu 
H 


Mode  2 


8 
10 
12 
14 
16 
18 
20 


•H 
CO 

e 

H 


Mode  3 


Terms    in  M 

4 

5 

6 

7 
11.258 

8 

9 

9.9897 

8.6493 

8.4420 

8.4173 

8.4157 

8 

.4153 

9.1333 

8.0673 

7.8967 

7.8763 

7.8748 

7 

.8746 

8.9990 

7.9743 

7.8060 

7.7857 

7.7840 

8.9927 

7.9683 

7.7987 

7.7777 

8.9750 

7.9537 

7.7847 

8.9347 

7.9237 

8.8897 

7.8913 

Terms   in  M 

4 

5 

6 

7 

8 

9 

9.1250 

8.5877 

8.5107 

8.5033 

8.5026 

8, 

,5026 

8.5133 

8.0630 

7.9977 

7.9913 

7.9906 

7, 

.9906 

8.4310 

7.9860 

7.9197 

7.9130 

7.9127 

8.4253 

7.9773 

7.9100 

7.9030 

8.4023 

7.9577 

7.8907 

8.3637 

7.9257 

8.3253 

7.8940 

Terms    in  M 

8 

a 

10 

c 

12 

•H 

14 

CO 

6 

16 

u 

CD 

18 

H 

20 

Mode  4 


15.592 

14.545 

14.371 

14.348 

14.347 

13.327 

12.616 

12.501 

12.487 

12.486 

16.288 

12.679 

12.046 

11.929 

15.959 

12.516 

11.897 

11.793 

15.910 

12.493 

11.872 

15.907 

12.488 

15.882 

12.470 

Terms 

;    in  M 

8 

16.418 

15.765 

15.660 

15.647 

15.646 

10 

13.714 

13.318 

13.258 

13.252 

13.252 

;s 

12 

14.761 

12.929 

12.589 

12.533 

c 

14 

14.463 

12.719 

12.388 

12.337 

co 

16 

14.415 

12.682 

12.350 

18 

14.412 

12.677 

it 
H 

20 

14.392 

12.662 
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the  ratio  of  hole  diameter  to  plate  width.   The  buckling  coefficients 
in  Table  1,  for  modes  1  and  2,  are  plotted  on  the  y  axis  (for  zero  d/b) 

The  values  for  modes  1  and  2  are  in  fair  agreement  with  experi- 
mental data.   No  correlation  between  the  theoretical  prediction  and 
the  experimental  observation  is  attempted  for  modes  3  and  4. 

The  buckled  deflection  of  the  plate  for  modes  1  through  4  is 
shown  in  Figures  5  through  8  respectively.   It  is  seen  in  these 
figures  that  the  four  possible  conditions  of  symmetry  are  satisfied 
by  the  four  modes.   This  result  is  given  in  Table  3. 

Table  3 
Resultant  Symmetry  for  the  Various  Mode  Shapes 

MODE  SYMMETRY  WITH  RESPECT      SYMMETRY  WITH  RESPECT 

TO  X  AXIS  TO  Y  AXIS 

1  f(y)  =  f(-y)  f(x)  =  f(-x) 

2  f(y)  =  f(-y)  f(x)  =  -f(-x) 

3  f(y)  =  -f(-y)  f(x)  =  f(-x) 

4  f(y)  =  -f(-y)  f(x)  =  -f(-x) 

c 
In  looking  at  Figure  4&,  it  is  apparent  that  mode  2  is 

converging  to  a  value  of  buckling  coefficient  that  is  somewhat 

higher  than  that  for  mode  1.   This  can  be  qualitatively  justified 

by  considering  the  symmetry  of  buckled  shape,  and  the  equivalent 

one  dimensional  problem. 

Consider  two  beams,  buckled  by  some  unspecified  loading,  into 

shapes  as  shown  in  Figure  9. 


Figure  9 
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It  is  known  from  beam  theory  that  the  shape  of  Figure  9b 
has  a  higher  critical  buckling  load  than  the  shape  of  Figure  9a. 
This  is  exactly  the  relationship  between  modes  1  and  2  of  the 
buckled  plate,  as  can  be  seen  in  Figures  5  and  6,  respectively. 

There  is  a  small  area  near  the  intersection  of  the  hole 
edge  and  the  y  axis,  which  assumes  cross  sectional  shapes  similar 
to  that  in  Figure  9a  for  mode  1,  and  in  Figure  9b  for  mode  2. 
It  is  only  in  these  small  areas  connecting  the  "top"  of  the  plate 
to  the  "bottom",  that  this  sort  of  bending  occurs,  so  that  the 
disparity  between  the  buckling  coefficients  of  the  first  and 
second  plate  modes  is  much  smaller  than  that  for  the  beam  modes. 

This  same  argument  applies  to  the  difference  between  modes 
3  and  4,  and  to  the  gross  differences  between  modes  3  and  4,  and 
odes  1  and  2. 


37. 


K 

8.5 

Q.k 
8.3 
8.2 
8.1 
8.0 
7.9 
7.8 
7.7 


k  $  6      7      8      9 

Number  of  terms  in  M 

FIGURE  k&   -  Plot  of  the  Buckling  Coefficient 
K  vs  the  number  of  terms  taken  in 
Mfor  mode  1. 
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8.6 


8.5 


8.1+ 


8.3 


8.2  1- 


8.1 


8.0 


7.9 


7.8 
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FIGURE  l|e  -  Plot  of  the  Buckling  Coefficient  K 
vs  the  number  of  terms  taken  in  N, 
for  mode  Ij.. 
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FIG   5  -   MAP   OF    DEFLECTION  FOR  MODE 
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FIG   6-  MAP  OF   DEFLECTION    FOR  MODE  2 
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APPENDIX  A 
SATISFACTION  OF  THE  BOUNDARY  CONDITIONS 
The  free  edge  boundary  conditions  for  the  hole  are  stated 
in  equations  15  and  16  as 


ww 
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The  plate  deflection,  W,  is  taken  as 
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where  the  coefficients  A>TW  and  B._,  are  varied  to  minimize  the  plate 

NM      NM  r 

energy,  and  the  coefficients  X^  and  Y   are  chosen  so  that  the 
boundary  conditions  expressed  by  equations  A-l  and  A-2  are  satisfied. 

The  derivatives  of  W  are  tabulated  in  Table  4.  When  these 
derivatives  are  substituted  into  equations  A-l  and  A-2,  and  evaluated 
at  r  =  1,  the  following  equations  result 
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where 
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S,  =  N(N  +  l)  -  N-0-  M2t) 

S-,  =  (N+  \)(N  +  Z)  -  (N+  I)  -0  -  M2t) 

S3  =  (N+2)(N+3)-(N+2)t)-mN) 


£  =  -N(N+|)(n+2)  +  N(N  +  i)+N  +  (2-^)M2N  +  (3-t>)m2 
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Equations  A- 4  and  A- 5  must  hold  for  all  values  of  9,  which  can 
only  be  true  if 

S,  +  SiXNM  +  S3YNM=  O  (A- 12) 


"£,+  CiXnm  +  't^YffM   ""  o 


(A- 13) 


for  any  value  of  N  and  M.   Thus,  a  value  of  X_„  and  Y._,  can  be 

NM      NM 

found  for  each  value  of  N  and  M,  such  that  the  boundary  conditions 
are  satisfied. 

By  application  of  Cramer's  Rule: 


y    = 

-S,     S3 
"t,      t3 

Anm 

S2    S3 

*a  ^3 

and 


'NM 


S2  ~S, 
tB  -t, 

ta  "ta 


*~  Sa"C<  "t"  S/"fc^ 
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The  above  results  are  incorporated  into  SUBROUTINE  BOUND  (see 
APPENDIX  F)  which  generates  values  of  X^  and  Y^  for  given  values 
of  N  and  M. 


52. 
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o 
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ii       u 
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I! 


J- 


u 

u        a> 

SB  > 


a> 
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APPENDIX  B 
EVALUATION  OF  THE  ENERGY  INTEGRALS 


The  total  energy  in  the  plate,  from  equations  25,  26,  and  27,  is 
Ut  =  U,+U2+U3 


where 


U      =  constant 


u2  = 


=  h. 


u3-| 


/-27T>oc 
'O     -'I 


[ovw  J  +  ^£  +  gjSfea&]  r  dr  da 


W~+^  +  Wg£ 


-2(|-T>) 


M^+W^-W 


(B-l) 


(B-2) 


drd9 


Note  that  these  integrals  contain  second  degree  terms  in  the 
derivatives  of  the  deflection,  W.   Since  W  is  represented  by  a 
series,  care  must  be  taken  to  insure  that  the  series  are  properly 
multiplied. 

The  product  of  two  double  sums,  both  summed  over  the  same  range 
for  the  subscripts  N  and  M,  is  a  quadruple  series,  with  dummy  subscripts 
I  and  K  summed  over  the  range  of  N,  and  K  and  L,  summed  over  the  range 
of  M. 

For  example; 


00  00 


OO  OO  OO    OO 


OP  OO 


OO  OO 


^^(i+rN)  cos  Me     =  Q^2](i+v-"I)cos  T9   V]T(i-hrK)cos  i_e 


1T-1  T=o 


K-t  L=0 


-V  Y  YY  (n-r_I)(i+r"K)cosxe  cos  Le 


1=1  J=0  K=l  L=0 


=VV y^(i+r"r+KK+r(r+K))cos  j©  cos  Le 
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It  is  seen  that  the  cross  product  terms  above  do  not  combine 

directly  into  one  term.   One  cross  product  term  involves  r  ^ the  other 

-K 
involves  r   .   For  this  reason,  cross  product  terms  in  the  energy 

integrals  B-l  and  B-2  must  be  identified  and  separated. 

As  an  aid  in  keeping  the  cross  product  terms  separated,  the 

notation  is  adopted  that  starred  (*)  functions  are  functions  of  I 

and  J,  and  that  non-starred  functions  are  functions  of  K  and  L,  in 

place  of  N  and  M,  respectively.   Equation  B-l  may  be  written 

-2TTV0O 


uz=  \\  J  [<rrvtf w.  +  ^£ 


+  •p:Trs(wr*W9  +  Wr%*)]  v-drdO 

and  with  the  substitution  of  the  Kirsch  solution  becomes 


(B-3) 


-£h_ 


-2TTV0O 


u,=^ 


'o   -/i 


-ia)wr*Wv+(l  +  !*-4)wJfWl.cos  2© 


r2 


r4    r* 


Equation  B-2   is  now 

U3=f/     /   (u3a  +  U3b+U3c)rdrd9 

where 


vd)-c/9 


u3b  =  d-o) 


w*(^+^w^%^ 


Wye  _  W6 

V  V-2- 


(B-4) 


(B-5) 


(B-6) 

(6-7) 
(B-8) 
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A   listing   of   the  expansion  of  the  products   of  the  derivatives 
necessary   for   the   evaluation  of  U_   and  U_   are   shown   in  Tables   5 
and  6,    respectively.      Further   presentation  of   the   energy   integrals    is 
in   the  notation  delineated   in  those   tables. 

Equations  B-4  and  B-5  now  become 

X27T 
[e1-E2+(ei+3E3-4E2]cos2©  +  E4.+Es  (b-9) 

-(E4  +  3^6)  cos  2©-(E70l4-E7b-3E9a-3E9b  +  2E9ft+EEgb)siM  ae]  de 


and 


,-— ,     r27T 


U3=  xXEZEj      [6.  +  62Q+62b  +  63a+63b  +  S5a+S5b  (B-1Q) 


+  S4.  +  S6  +  (I-t))  6ea+S2b+S3a+  63b 


+  2(i-t))  S7  +  69  +Gea  +  G8b 


d0 


When   the   integration  in  ©  is   done,         equations  B-9  and  B-10  reduce   to 


oo   oo    oo   oo 


"An- 

AKL" 

0(r 

_BXJ 

B«L 

(B-ll) 


and 


OO      QO      OO      OO 


r-i  3=0  k-i  \.=o 


Au-AKl 

ov- 
Bij'Bkl 


(B-12) 
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Tables  5  and  6 

Tabulation  of  the  Integrals  of  the  Radius  Dependent  Terms 

NOTATION: 

G1 ,  G9 , . . .E  ,  E  , ....  etc.  are  functions  of  I,  J,  K,  and  L. 


means  5L-.  as  used  previously,  only  as  a  function  of  K  &  L. 
x*  means  X   as  used  previously,  only  as  a  function  of  I  and  J, 


x 


y  means  Y,„,  as  used  previously,  only  as  a  function  of  K  and  L. 

y  means  Y,„,  as  used  previously,  only  as  a  function  of  I  and  J. 

J  NM         r 


THE  INTEGRALS  ARE  TABULATED  ON  THE  FOLLOWING  PAGES.   THE  FIRST 
NINE  TERMS  IN  EACH  COLUMN  ARE  SUMMED  TO  FORM  THE  WHOLE  EXPRESSION, 
THE  TERMS  a(J,L)  or  b(J,L)  SHOULD  BE  CHOSEN  FOR  THE  UNKNOWN 
COEFFICIENTS  A^  or  B   RESPECTIVELY. 
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Table  £ 

Tabulation  of  the   Compacted  Farm  of  the  Terms   in  the 
Membrane  Energy  Integral 


h  = 


E, 


E. 


\ 


Er 


E7. 


r  r 


w  w 

r  r 


r  r 


3T 


.00 


wew© 


CO 


W   W" 


^b 


©  r 


oo 
©  r 


i 


dr 


dr 


dr 


dr 


dr 


dr 


dr 


=  f . 


Lb(J'L)'BIJ,BKLJ 


=  f , 


Lb 


=  f . 


=  f . 


=  f , 


L* 


=  f 


6'i 


=  f 


7a 


a 
b 


f7b|. 


(J,L) 

or 

(J,L) 

(J,L) 

or 

(J,L) 

(J,L) 

or 

(J,L) 

(J,L) 

or 

(J,L) 

(J,L) 

or 

(J,L) 

(J,L) 

or 

(J,L) 

(J,L) 

or 

( J,L) 

(J,L) 

or 

(J,L) 

lij,akl] 


A 


IJ>AKL 


BIJ,BKLj 


■  \l 


BIJ 


U 


B 


IJ 


A 


B 


IJ 


IJ 


A 


IJ 


B 


IJ 


LIJ 


B 


IJ 


V 
BKlJ 


TSLi 

5klJ 


BKL. 

akl] 
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Table   5>    (continued) 


3e. 


Co 


K  - 


dr 


=  f 


[aU.U-Ajj'Aja 


8a 


or 


bCJ.U.Bjj-B^j 


E8b  " 


00 
9  r 


dr 


=  f 


8b 


aU.D-Ajj-AjQ.-! 

or  j 


E 


9a 


e  r 


dr 


=  f 


;a(J»L>-Au'AKL 


9a! 


or 


jbtJ.D.Bjj-B^j 


'9b 


OO 


e  r 
"3" 


dr 


=  f 


9b* 


^(^D.Ajj.A^-j 


or 


lb(J,L)  -Ajj-AqJ 
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Table  5  (continued) 


*»        -  IK 


+      KK+l)x 
(I+K+l) 


IlK+2)y 

(I+K+2) 


+      (I+l)Kx 

(I+K+l) 


+ 


+ 


(I+l)(K+l)xx 
(I+K+2) 

(I+l)(K+2  )jx 


(I+K+3) 


■f         — 

IK 

r2 

(I+K+2) 

+ 

I(K+l)x 
(I+K+3) 

+ 

I(K+2)y 

+ 

(l+l)Kx* 
(I+K+3) 

4. 

(I+l)(K+l)xx* 

(I+K+ij.) 

+ 

(I+l)(K+2)yx* 

+      (l+2)Ky'"'  (I+2)Ky" 

(I+K+2)  (I+K+I4J 


where  where 

a(J,L)   =  cos    J©   cos   L©  a(J,L)   =  cos   J©  cos   L© 

b(J,L)   =  sin   J©  sin  L©  b(J,L)   =  sin  J©  sin  L© 
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Tabl©  5  (continued) 


IK 


(I+K+lj.) 
I(K+l)x 


KK+2) 

TT+K+5 


r 


(l+l)Kx 


(l+l)(K+l)xx'' 
+     (I+K+6) — 

+  (l+K+?) 


* 


(3>2)Ky'"" 
(l+K+61 


+   (I+2)(K+l)xy'f 
(I+K+7) 

+   d+2)(K+2)yy'"" 

+   (i+k+8) 


JL 


+ 


TT+KT 


JLx" 
+  (i+K+i) 


JLy 


(I+K+2) 


(I+K+l) 


,    JLxx" 

+    (i+k:+2) 


4.     JLxy'"" 

(I+K+2) 

(I+k+3) 

(I+K+lj.) 


where 

a(J,L)  =  cos  J9  cos  L9 
b(J,L)  =  sin  J9  sin  L9 


where 

a(J,L)  =  sin  J©  sin  L9 

b(J,L)  =  cos  J9  cos  L9 
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Table  £  (continued 


H  m 

JL 
(I+K+2) 

+ 

JLx* 
(I+K+3) 

+ 

JLy* 
(I+H) 

+ 

JLx 

(I+K+3) 

+ 

JLxx" 
(I+K-Hj.) 

+ 

JLxy" 
( I+K+5 ) 

+ 

(T+l+U 

+ 

JLyx'"" 

+ 

[iMh\ 

where 

a(J,L) 

=  sin   J9   sin  L9 

b(J,L) 

=  cos    J9   cos   L9 

+ 


JL 


I+K+i| 
JLx* 

JLy'"' 


I+K+6 


JLx 


T+it+F 
JLxx" 


I+K+6 


JLxy 


I+K+7 

JLy 
I+K+6 

JLyx" 


I+K+7 

JLyy" 
I+T+S 


where 

a(J,L)  =  sin  J9  sin  L9 

b(J,L)  =  cos  J9  cos  L9 
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Table  5  (continued) 


7a 


"  TT+kT 

J(K+l)x 
(1+K+iJ 


f  7b  : 


LI 

Ti+kT 


L(I+l)x~"' 
( I+K+l) 


J(K+2) 
( I+K+2 


f 


L(I+2)yw 
(I+K+2) 


JKx'"' 
( I+K+l J 

J(K+l)xx" 
(I+K+2) 


J(K+2)yx""" 
(I+K+3) 

JKy  " 
(I+K+2) 


J(K+pxy* 
"    (I+K+3) 

J(K+2)yy* 
"    (I+K+l[) 

where 

a(J,L)  =  -sin  J9  cos  L9 

b(J, L)  =  sin  L9  cos  J9 


LIx 
(I+K+l) 

L(I+l)xx': 
'"(I+K+2) 


L(I+2)xy'"' 
-    ( I+K+31) 

Lly 

(I+K+2) 

L(I+l)yx* 
(I+K+3) 

.  L(I+2)yy"' 
(I+K+l;) 

where 

a(J,L)  =  -sin  L9  cos  J9 

b(J,L)  =  sin  J9  cos  L9 
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Table  5  (continued) 


8a 


JK 

TT+K+2 ) 


(i+K+3) 

J(K+2)y 
(I+K+I4.) 


8b 


LI 
'    (l+K+2) 

L(I+px" 
(I+K+3) 


L(I+2)y'' 


JKx"" 
(I+K+3) 

J(K+l)xx' 
( I+K+ij.) 


J(  K+2  to- 
ll+K+5) 

JKy" 

rr+fen 

J(K+l)xy* 
Cl+K+5) 

J(K+2)yy* 
(I+K+oV 


where 

a(J,L)   =  -sin  J9   cos   L© 

b(J,L)   =  sin  L©  cos    J© 


Llx 
(I+K+3) 

L(I+l)xx'::" 
(I+K+ij.) 

.     L(I+2)xy* 
Tl+K+J) 

Lly 
(I+K+ij.) 


L(I+l)yx'" 
(I+K+5) 


.     L(I+2)yy" 
(I+K+6) 

where 

a(J,L)   =  -   sin  L©  cos    J© 

b(J,L)   =  sin  J©   cos   L© 
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Table  5  (continued) 


9a 


JK 
(I+K+Ij.) 


J(K+l)x 

(I+K+5) 


J(K+2)y 
(I+K+6) 


JKx'"~ 

rx+K+n 


J(K+l)xx""" 
(J+K+6) — 


J(K+2)yx'"' 

(I+K+7) 


JKy'"' 
"   (I+K+6) 

J(K+l)xy* 

J(K+2)yy"" 
'   (I+K+8V 

where 

a(J,L)  =  -sin  J9  cos  LO 

b(J,L)  =  sin  L9  cos  J9 


9b 


LI 

"  (I+K+Ij.) 

L(I+l)x* 

(I+K+5) 

L(I+2)y* 
(1+K+ol 


LIx 

TT+1T+TT 

L(I+l)xx': 
(I+K+6) 


L(I+2)xy 
( I+K+7 5 


Lly 
(I+K+6) 


L(I+l)yx" 
(I+K+7) 


(I+K+8) 


where 

a(J,L)  =  -sin  L9  cos  J9 

b(J,L)  =  sin  J9  cos  L9 
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Table  6 

Tabulation  of  the  Compacted  Form  of  the  Terms  in 
the  Bending  Energy  Integral 


r-OO 


W  VT  dr 
rr  rr 


2a 


oo, 


r  rr 


dr 


-  o 


2b 


oo 

W  W* 
r  rr 


;    r 
'  o 


dr 


co 


3a 


3b 


a, 


WrrW99 


oo 
/   rr  9© 


Jo       r 


co 

C  W  W 

r  r 


dr 


dr 


dr 


'$* 


oo 

W99Wr 


dr 


-', 


G 


5b 


oo 

r  w   w 


dr 


"  Si 


:<a 


=  g 


2a 


:b 


=  g 


2b 


=  g 


3a 


'3b 


a 

I? 


=  gi 


H 


g5a 


L* 


S5b 


(J,L] 
or 

'AIJ' 

(J,L] 

'BIJ' 

(J,L] 
or 

*AIJ* 

(J,L] 

'BIJ* 

(J,L] 
or 

[*AIJ' 

(j,l: 

>'BIj' 

or 

|,AIJ' 

(j,l: 

,,BIJ' 

or 

,<AIJ' 

(j,l: 

''BIJ' 

(J,L 
or 

>,AIJ' 

(J,L 

>'BIj' 

(J,L 
or 

>*AIJ* 

(J,L 

>'BIJ' 

(J,L 
or 

>'AIJ' 

(J,L 

^IJ* 

hL   j 

AKL 
BKL. 


B. 


KL 


AKL 


B, 


KL 


akl] 

BKLJ 


BKLJ 


akl] 


66 


Table  6  (continued) 


oo 


weowee 


3T 


dr 


(b(J,L).B1J.BKL 


OO 


W     ¥" 
r0  r© 


dr 


=     g 


7  ' 


"a(^;L)"Aij"AKL" 

or 


oo 


G 


8a 


W0WrQ 


,3 


o 


dr 


;a(  J,L)-ATT-A-T-] 

S8a'j      or  I 


oo 


r8b 


i    w  w" 
/      ©  r© 


J       * 


-  dr 


ja(J,L).A      -A^" 
go-u.-        or 


OO 


©W© 


1 


3T 


dr 


=     S< 


aUD-A^-A^ 

or 
b(J,L)-BIJ.Bn.j 
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Table  6    (continued) 


R       =     IK(I+1)(K+1) 
sl  (l+K+2) 

+      (I+D (K+l) (1+2) (K+2)xx* 
(I+K+lj.) 

+      (I+2)(K+2Hl+3)(K+3)yy* 

(i+K+6) 

.      K(K+l)(I+l)(I+2)x* 
(I+K+3) 

.      I(I+l)(K+l)(K+2)x 

+    hc+K+TJ 

.      K(K+l)(I+2)(I+3)y* 

(i+K+4) 

+     Id+1)  (K+2)(K+3)y 
(I+K+4) 

+      (I+l)(I->-2)(K+2)(K+3)yx^ 
(I+K+5) 

,      (K+l)(K+2)(I+2)(I+3)xy* 

+    (i+K+5) — 

where 

a(J,L)   =  cos    J©  cos   L© 

b(J,L)   =  sin  J©  sin  L© 
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Table  6    (continued) 


g 


IK(K-H) 
2a  '     ~    (I+K+2) 


I(K+l)(K+2)x 
(I+K+3) 

I(K+2)(K+3)y 
(I+K+lj.) 


g 


2b 


IK(I+1 
(I+K+2 


K(I+l)(I+2)x 
(I+K+3) 

K(i+2)(l+3)y 

(l+K+C) 


tt 


(I+l)(K+l)Kx'"' 
(I+K+3) 

(I+l)(K+l)(K+2)xx* 

TT+TRn 

(I+l)(K+2)(K+3)yx' 

— rr+K+in 


(I+2)(K+l)Ky' 
"      U+K+lj.) 


(I+2)(K+l)(K+2)xy' 
(I+K+5) 

(I+2)(K+2)(K+3)yy' 
(I+K+6) 


(K+l) 1(1+1 )x 
(I+K+3) 

(K+l)(I+l)(I+2)xx* 
(I+K+U 

(K+l) (1+2) (1+3 )xy* 

— (i+KA) 

(K+2)I(I+l)y 
(I+K+I4.) 

(K+2)(I+l)(I+2)yx*' 
(I+K+5) 

(K+2)(I+2)(I+3)yy* 
(I+K+6) 


where 

a(J,L)   =  cos    J©  cos   L© 

b(J,L)   =  sin  J©  sin  L© 


where 

a(J,L)   =  cos   J9  cos   L© 

b(J,L)   =  sin   J0  sin  L© 
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Table  6    (continued) 


g 


3a 


(I+l)(I+2)L2x* 


I(I+1)L2 
'      (I+K+2) 

I(I+l)L2x 
(I+K+3J 

I(I+l)L2y 

l)(I+2 
(I+K+3 


(I+l)(I+2)L2xx* 
(I+K-H4J 

(I+l)(I+2)L2yx* 

TTSB5T 

(I+2)(I+3)L2y* 
(I+K-Hj.) 


where 

a(J,L)   =  cos    J9  cos   LO 

b(J,L)   =  sin  J6  sin  LO 


g 


3b 


KCK+pJ2 
(I+K+2) 


-        \i+K+3) 


K(K+1)J 


2* 


(I+K+1+1 
(K+l)(K+2)J2x 

-     (i+g+ir — 


-:c- 


„     (K+^HK+gj^yx* 

.     !K+2l(K+3)jfy£_ 
(I+K+6) 

where 

a(J,L)   =  cos   J©  cos  L9 

b(J,L)   =  sin  J9  sin  L© 
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Table  6  (continued) 


6. 


h 


IK 

n+ir+2) 


.    K(in)x'' 
+     (I+3K+3) 


K(I+2) 
(I+K+J 


a- 


& 


.       (K+l)Ix 

+     (l+K+3) 

(K+l)(I+l)xx" 
+         (liK+t) 

,      (K+l)(I+2)xy 
(I+K+5) 

+      jK+2)Iy 


,       (K+2)(I+l)yx, 
+     — (l+K+5) 


+     (K+2)(I+2)yy': 
+     (I+K+6)  ^ 


SS 


a 


j2L 

TTTkT^t 

J^K+l): 

jfOM)y 
(I+K-UjT 

Ax* 

(1*3) 


(I+K+i) 


J2(K+2) 


ff 


(I+K+ 


J*Ky* 

(l+k+5) 


«r(K+2) 

(I+K+6 


f~ 


where 

a(J,L)   =  cos   JO   cos   L© 

b(J,L)   =  sin   J©  sin  LO 


where 

a(J,L)   =  cos   J©  cos   L© 

b(J,L)   =  sin   J9   sin  L© 
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Table  6    (continued) 


_        L2I 

*6     = 

J2!2 
(I+K+2) 

.      L2(I+l)x* 
+        < I+K+3) 

+ 

Aftt* 

( I+K+3) 

.      L2(I+2)j* 

+      (T+M1 

+ 

d+A) 

.      L2Ix 
(I+K+3) 

+ 

d+K+3) 

L2(I+l)xx'"" 
(±+K+i) 

+ 

(I+K+i;) 

TT+k+5) 

+ 

J2L2xy* 
(I+K+£) 

+     (I+^A) 

+ 

j2L2y 
(I+A) 

.      L2(I+l)yx* 

+ 

J^yx* 
(I+K+5) 

,      L2(I+2)yy* 

+ 

jVyf 
(I+K+o) 

where 

where 

a (J, L)   =  cos    J9  cos   L© 

a(J,L)   =  cos    J©  cos   L© 

b(J,L)   =  sin  J©  sin  L© 

b(J,L)   =  sin  J©  sin  L© 
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Table  6  (continued) 


g- 


IJKL 
(I+K+2) 


IJ(K+l)Lx 
(I+K+3) 


«8i 


JK 
(I+K+2) 

J(K+l)x 
(f+K+3) 


.   IJ(K+2)Ly 
+   (1+K+C) 

(I+l)JXLx'"~ 

+  (i+k+3) 

(I+l)J(K+l)Lxx" 
CT+K+I4.) 

,   (I+l)J(K+2)Lyx* 
+  (i+K+S) 


■$:- 


(1+2) JKLy 
(I+K+tj.) 

(I+2)J(K+l)Lxy 
(I+K+5) 

(l+2)J(K+2)Lyy' 
(I+K+6) 


* 


J(K+2)y 
Tr+K+IT) 


JKx" 
(I+K+3) 

J(K+l)xx*' 
(I+K+I+) 


J(K+2)yx"' 
(I+K+5) 


(I+K+i^) 
J(K+l)xy* 


J(K+2)yy 
(I+K+6) 


# 


where 

a(J,L)  =  sin  J©  sin  L© 

b(J, L)  =  cos  J©  cos  L© 


where 

a(J,L)  =  sin  J©  sin  L© 

b(J, L)  =  cos  J©  cos  L© 


73 


Table  6  (continued) 


*8i 


LI 
"  (I+K+2) 

L(H-l)x* 
(l+K+3) 


g 


■(1+22* 

(I+K+I+ 


# 


■3:- 

r 

LIx 
(I+K+3) 

L(I+l)xx 

Ct+rt+3)' 


L(I+2)xy" 
Lly 

-  TirrajT 

L(I+l)yx* 
"   (I+K+51 

L(I+2) 


T 


I+K+ol 


where 

a(J,L)  =  sin  JO  sin  L9 

b(J,L)  =  cos  J0  cos  LO 


:9 

JL 
(l+rf+2) 

+ 

JLx 
(I+K+3) 

+ 

JLy 
(I+K+lj.) 

+ 

JLx 

(-1*4) 

+ 

JLxx 

(MM*)' 

+ 

JLxy" 
(I+K+5) 

+ 

JLy 

+ 

-rr    * 

JLyx 

(T+K+5) 

+ 

JLyy" 

(I+K-rfe) 

where 

a(J,L)  =  sin  JO  sin  LO 

b(J,L)  =  cos  J©  cos  L© 
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APPENDIX  C 
INTEGRATION  OF  SINE-COSINE  TRIPLE  PRODUCTS 

When  evaluating  the  energy  integrals  for  the  plate,  it  is 
necessary  to  evaluate  the  following  two  integrals 
-27T 


C 
Ic=  /   (cos  A9  cos  BO  cos  C&U& 


and 


•27T 


T    -  I        sin  A©    sin  B©  cos  C©    d© 


(C-l) 


(C-2) 


where  A,  B,  and  C  are  integers. 

The  above  integrals  can  best  be  evaluated  by  expanding  the 
integrands  using  the  identities 

Sin  y.   cos  y  =  ■£  siN(-x+y)+  SiN(**-y)| 


COS  7.  SIN  y  -  L fsiM(x+t|)-SiN('X.-y)j 
cos  x  cosy=.  ^-[cos(*+ij)-»-co$(,x-y)] 
Sin  x  sin  y  =  -i-fCosU-y)-  c°s("*+y)l 

The  object  is  to  reduce  the  integrals  to  the  form 


(C-3) 
(C-4) 
(C-5) 
(C-6) 


-2TT 

sin  D©   d© 


or 


I 


znr 

cos   D©  d© 


(C-7) 


-o  -o 

where  D  is  an  integer  sum  or  difference  of  A,  B,  and  C.   Note  that  the 
sine  integral   in  C-7  is  always  zero  and  that  the  cosine  integral  is 
nonzero  only  if  D  is  zero. 

Expanding  equation  C-2,  using  identity  C-3,  we  get 


rs=  £/  sm  a© 


siN(e+-c)©->-  sin(b-c)9 


d© 


(C-8) 
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Continuing  the  process  with  identity  C-6 
I$  =  4/  fcos(A-B-C)ft-cos(A+B+c)e4-cos(A-fc+c)e-cos(A*-B-c)© 


de 


(C-9) 


I   is  now  expressed  as  the  sum  of  four  integrals 


Is=I|  +  l2+l3  +  U 


where 


(C-10) 


Xl  = 


h   = 


If/2      if  (A-B-C)  =  0 
0  otherwise 

-IT/2   if  (A+B-fC)  =  0 
0  otherwise 


I3  = 


h  = 


If/2      if  (A-B-K3) 

0  otherwise 

-If/ 2      if  (A+B-C) 

0  otherwise 


=  0 


=  0 


A  similar  expansion  for  I   (equation  C-l)  can  be  made.   The 
results  are  the  same  as  above,  except  that  the  signs  of  I?  and  I, 
are  positive. 
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APPENDIX  D 
SOLUTION  OF  THE  EIGENVALUE  PROBLEM  (Rx  +  T)  A  =  0 

The  canonical  form  of  the  matrix  eigenvalue  problem  is 

(iX-g)X^O  (D-D 

where  I   is  the  identity  matrix,  X  an  unknown  eigenvector,   \  ia 
an  unknown  scalar  parameter  (eigenvalue),  and  B  a  known  coefficient 
matrix. 

Equation  D-l  has  a  solution  if  determinant  of 

(IX-|)  =  0  (D-2) 

If  this  determinant  is  expanded,  a  polynomial  equation  in  K   results. 
However,  if  the  determinant  is  of  order  greater  than  three  or  four, 
it  is  generally  inefficient  to  obtain  this  polynomial  equation  and 
then  solve  it. 

There  are  alternative  procedures  for  solution  of  the  eigenvalue 
problem.   A  method  which  gives  extremely  accurate  eigenvectors  and 
eigenvalues  for  symmetric  matrices  is  one  named  for  Jacobi,  which 
is  presented  by  Fox  (1).   Brown  and  White  (13)  give  a  comparison  of 
the  accuracy  of  the  Jacobi  method  with  that  of  other  methods  available, 
The  necessity  for  extremely  accurate  eigenvectors  is  shown  later. 

Consider  the  eigenvalue  problem  (equation  33) 

(EX*-X)A=o  (D"3) 

where  both  R  and  T  are  symmetric  matrices  and  T  is  positive  definite. 
Dividing  equation  D-3  by  X  ,  and  making  the  substitution 

3B-  J  (D-4) 

gives 
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Any  symmetric  matrix  with  real  elements  can  be  reduced  to  a 
product  of  the  eigenvector  and  eigenvalue  matrices,  so  that 

X  =  iRiT  (D-6) 

where  D  is  a  diagonal  matrix  of  the  eigenvalues  (d.)  of  T,  and  Q 
is  a  matrix  of  normalized  eigenvectors  (q.),  obtained  from  a 
solution  of  the  problem 

Note  that  ^  =  £,  since  the  matrix  T  is  symmetric  and  has  orthogonal 
eigenvectors . 

Equation  D-5  may  now  be  written 

(§£gT<j  ~E)A=  0  (D-8) 

Premultiplying  by  Q~  ,  we  get 

(2§T^-girR)A  =  o  (D-9) 

T 
Factoring  a  Q  out  of  A  =  JA  =  ^  ^  A,  and  letting 

Y  =  QTA  (D-10) 

the  problem  reduces  to 

(gy-gTg§)Y  =  o  (D-ii) 

Since  T  is  positive  definite,  D  contains  all  positive  terms. 
By  again  premultiplying. this  time  by  D   2,  factoring  D  2,  and  making 
the  substitution 

\m^l  (D-12) 

equation  D-ll  reduces  to 
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and  becomes 

which  is  of  canonical  form. 

%     -h 
Note  that  D  and  D   are  easily  obtained  by  taking  the  square 

roots,  and  the  square  root  reciprocals,  respectively,  of  the  diagonal 

terms.   Also  note  that  15  is  symmetric.   The  problem  is  now  well 

suited  for  solution  by  use  of  the  Jacobi  method.   It  is  seen  that 

inaccuracies  in  obtaining  eigenvectors  Q   would  cause  compounded 

errors  in  the  final  solution  for  A  and  X   • 

The  problem  is  returned  to  original  form  with  the  substitutions 

X  =  ~^  (D-15) 


and 


A  =  QD"'/aZ  (D-16) 
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APPENDIX  E 
SUGGESTIONS  FOR  EXTENSION  OF  THE  PROBLEM 

There  are  two  obvious  directions  for  extension  of  this  problem 
of  a  thin  plate  with  a  circular  hole.   The  buckling  problem  for 
a  strip  of  finite  width  and  the  postbuckling  problem  for  both  finite 
and  infinite  plates  remain  to  be  treated  analytically. 

The  postbuckling  problem  is  the  more  difficult  of  the  two 
because  of  the  coupling  between  the  stress  distribution  and  the 
buckled  deflection.   Johns  (2)  outlines  the  postbuckling  situation 
for  plates  and  suggests  several  methods  of  solution.   Solution  of 
the  postbuckling  problem  would  be  important  because  of  the  revised 
stress  concentration  factors  that  would  be  obtained;  cf.  Houghton 
and  Rothwell  (16) . 

Extension  of  this  problem  to  a  finite  strip,  of  width  2b,  would 
require  the  following: 

a.  Using  the  Howland  solution  (17)  instead  of  that  due  to  Kirsch, 
to  obtain  the  stress  distribution  in  the  plate. 

b.  Satisfying  additional  boundary  conditions  for  two  "free" 
straight  edges. 

Integration  of  the  energy  terms  for  the  finite  strip  is  between 

the  limits  r  =  <X   and  r  =   — ,  so  that  the  integration  with 

SIN   Q 

respect  to  9  becomes  more  complicated  than  for  the  infinite  plate. 

This  integration  is  over  the  interval  0  to  27T.  Wylie  (10) 
presents  a  method  for  the  evaluation  of  real  trigonometric  integrals 
over  this  interval.   The  method  involves  substitution  of  appropriate 
complex  variables  for  the  sine  and  cosine  terms,  and  then  evaluation 
of  the  residues  inside  a  unit  circle. 
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There  are  two  possible  methods  for  satisfaction  of  the 
additional  boundary  conditions  of  the  finite  strip.   Introduction 
of  additional  unknowns,  similar  to  X^  and  Y   in  equation  A-3, 
into  the  assumed  series  for  deflection  could  be  made.   These 
unknowns  could  be  determined  by  simultaneous  solution  of  all  the 
boundary  condition  equations.  Another  method,  the  Lagrangian 
multiplier  method,  introduces  dummy  variables  so  that  additional 
restraining  equations,  those  expressing  the  boundary  conditions, 
can  be  solved  simultaneously  with  the  energy  expressions. 
Budiansky  and  Hu  (12)  present  a  clear  example  of  the  use  of  this 
method  in  minimizing  the  plate  energy,  while  satisfying  plate 
boundary  conditions. 
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APPENDIX  F 
COMPUTER  PROGRAM  LISTINGS 

The  following  pages  contain  a  listing  of  the  computer  programs 
used  in  the  solution  of  the  problem.  All  programs  are  written  in 
FORTRAN  IV  for  use  on  an  IBM  System  360  computer. 

The  control  card  listing  on  the  first  page  is  effective  as  of 
September  1967.   Users  of  this  program  should  check  current  procedures 

The  maximum  terms  to  be  taken  in  N  and  M,  and  the  plate  thickness 
•  read  into  the  computer  by  data  card  (see  comment  cards  regarding 
input  requirements).   Young's  modulus,  E,  and  Poisson's  ratio  are 
established  by  cards  serial  0010  and  0011  in  the  main  program. 

Printed  output  from  the  computer  is  a  column  matrix  EIVU  which 
is  an  array  of  all  of  the  values  of  S   in  psi  for  the  parameters 
of  plate  thickness,  Poisson's  ratio,  Young's  modulus,  terms  in  N, 
and  terms  in  M.  Matrix  EIVU  is  followed  by  a  map  of  the  buckled 
plate  deflection  for  the  minimum  value  of  tensile  critical  buckling 
stress . 
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S   =  KE(t/a)2 
cr 

where  E  is  Young's  modulus  of  elasticity,  t  is  plate  thickness,  a  is  a  hole  radius, 

and  K  is  a  buckling  coefficient  having  the  following  values  for  the  first  four 

modes:   7.65,  7.80,  11.7,  12.3.   A  comparison  is  made  with  some  experimental 

determinations  of  Danis  for  plates  of  finite  width. 
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